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It is shown that the field strength according to mode theory and ray theory in the 
VLF band are derivable from the same expression of the original vector potential, and the 
result of one theory is the analytic continuation of the other one in another range of con- 
vergence. In fact, both ranges of convergence overlap. Estimates of these ranges are 
made and an example shows that within this overlapping region (between distances of 300 
and 2,000 km) both theories give the same result. Using this fact calculations of frequency 
spectra are possible which in the case of a white noise dipole show some similar features 
to measured frequency spectra of lightning discharges. 



1. Introduction 

Electromagnetic waves at frequencies below 100 
kc/s are reflected during daytime at the bottom of 
the ionosphere at about 70 to 80 km lieight. At 
this height the collision number is still greater than 
the angular frequency co, and the depth of penetra- 
tion of these waves into the ionosphere is limited. 
For distances greater than about 500 km and fre- 
quencies below about 50 kc/s tlie actual ionospliere 
during the day can be replaced to a good order 
of approximation by a homogeneous isotropic plasma 

with conductivity 0-^= — ^— and a sharp lower 

boundary. {N ^ and v are mean values of electron 
density and collision number respectively, and at is 
of the order 10"^ 12"^ cm~^) 

For the ''medium'' distaiu'cs between about 500 
km and the limiting geometric skip distance of the 
first hop wave (1,700 km at 70 km reflection height) 
higher order terms both in mode theory and ray 
theory may be important. It is the purpose of this 
paper to compare mode theory and ray theory in 
this distance range. 

2. Setting up the Equations 

We wish to calculate the vector potential of the 
commonly assumed vertical dipole within the system: 
flat earth and isotropic homogeneous ionosphere 
with sharp lower boundary. Following the method 
of H. Weyl [1919] it can be shown that the vector 
potential of this vertical dipole is 

-> 



1 Contribution from Heinrich-nertzlnstitut, Berlin-Charlottenburg, Germany 
(Western Zone). 



+ A'f6-^'^^^2n^+^o+^^ \1 0^3^A (1) 



with 



Bi- 



'n\C+Jn\-S^ 

e \ e 



ni^l- 



ecco 



osphere 
th 



(7, conductivity of -< 
e \ car 

(p,4>,z) cylindrical coordinates 



} height of the antenna of-< ^ 
^ \^ transmitter 

h virtual height of the ionosphere 



p=Po e^"' moment of the dipole 
/:=— wave number 

A 

The flrst two terms of the integral (1) represent the 
vector potential of a vertical dipole above a con- 
ductive earth without taking into account the 
ionosphere, a well-known problem that has been 
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solved by A. Sommerfeld [1948] and which gives 






(2) 



tF(Q:) = |Pr|e ^^^ is the Sommerfeld absorption factor 
a——jkp[l—ne{l+nl)~'^^^ is the numerical distance. 

Again, following Weyl [1919], the terms of the sum 
within the integral are approximately 
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This (saddle-point) approximation is valid when 
Vyn^X and Cyn is not too small compared with unity. 
\ni\ should also be somewhat greater than unity. 

Setting z=Zo=0 and neglecting higher order 
terms of Tyn the vertical component of the electric 
field strength becomes 
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=2E„[ 



W{a) + 2^sm'd„ 



RfidjRr'ien) 



{i+RAe„))\ 



jfc(p-r 



->] (4) 



with £^0= "" (free space value of field strength 

of the vertical dipole at the distance p) 



sin dn = -- 
' n 



rl=p^+{2nhy 

Equation (4) was derived explicitly by Wait and 
Murphy [1957], using a geometrical-optical approach. 
It can be interpreted in the language of ray-optics 
as consisting of the ground wave, a wave reflected 
from the ionosphere, and further waves reflected 
between the groimd and the ionosphere any number 
of times. 

Carrying out the integration over 0' and forming 
the sum, (1) becomes 
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Following Budden [1951] and Alpert [1955] (who 
set Re=l) equation (5) may be approximated by a 
sum of modes. For s=Zo=0, it has the form 



n(p,</.,o)= 
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Here the values Cn are determined from the eigen- 
value equation 



i^^(OJg-2^■''^^« = 6-2^■^^ 71 = 0, 1, 2, 



(7) 



The vertical component of the electric field is then 
(using the approximation of the Hankel function 
H^^\x) valid for large values of x) 
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for />5 kc/s and —^ is of the order of unity. 

Equation (8) in the language of wave optics con- 
sists of the possible propagation modes. Actually it is 
the approximate continuation of the ray equation into 
another range of convergence. More general forms 
of eqs (6), (7), and (8) are given by Wait [1957a] 
and Wait and Spies [1960] for an imperfectly con- 
ducting and curved earth. 

3. Estimate of the Ranges of Convergence 

We wish to confine the sums in eqs (4) and (8) to 
the first five terms. This is reasonable provided 
that the sixth term is smaller than 1 percent of the 

first one. Thus, for— ^=1, we get from eq (4) 

eoco 
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or 
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From eq (8) we find 
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Equation (11) is the first approximation of the 
eigenvalue-equation (7) according to J. R. Wait 
[1957a]. 

Because we have 



«1 



eq (10) becomes 
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The range of small convergence of eq (8) simul- 
taneously is the range of good convergence of eq (4) 
and vice versa. In fact, both ranges overlap. 

4. Comparison Between Ray Theory and 
Mode Theory 

Equations (4) and (8) give only approximate 
results of equation (1). To get an idea of their 
accuracy examples have been calculated with the 
following parameters: 



i?. = l 



eoco 



:1; y=15kc/s. 



For the determination of the eigenvalues Cn. (eq (7)) 
tables have been used which are reported elsewhere 
[VoUand, 1959]. Table 1 contahis values of the 



argument 



and the phase of 



2^0 



depending on the distance p with an ionosphere of 
constant height h=70 km, calculated from eq (4) 
and eq (8), respectively. Table 2 shows the same 
values depending on the lieight h and with constant 
distance p= 1,000 km. Figure 1 is a graphical form 
of tables 1 and 2. 

From tables 1 and 2 it can be seen that between 
300 and 2,000 km distance and between 35 and 100 
km height, respectively, both equations give the 
same result to within about 5 percent in argument. 
Figure 2 illustrates the influence of the higher order 
terms for the case h = 70 km. In figure 2a, corre- 
sponding to eq (4), the distance 0-A represents the 
constant ground wave. The vector drawn from A 
to one of the circles represents the first liop wave, 
the thhi lines being equivalent to the vectors of 
second, third, and fourth hop waves, respcH-lively. 
The vector drawn from the point O to the solid line 
which connects the end points of the thin hues gives 
the values of argument and phase for the different 
distances tabulated in table 1, second and fourth 
row. 

Figure 2b shows the first five terms of the sum 



n=l \Si/ 



Jkp{S,-S^) 



in the same sequence as in figure 2a. Tlie distance 
0-A is equal to the constant first term 1, etc. The 
vector (h'awn from the point O to the solid line must 
be inulti|)lie(l hy the first term of eq (8) 



g-,V/4^3/2g,-fcp(l-S0 



_pX 

h 



to get the values in the first and tliird r'ow of table 1. 



Table 1. Argument a and phase p of Ez/2Eo as a function of the distance p with constant h=70 km according to equation (4) (symbol 

'V") and equation (8) (symbol "w^'). 
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Table 2. Argument a and phase p of EJ,2Eq as a function of the height h with constant p = 1,000 km according to equation {4) 

(symbol 'V") and equation (8) (symbol "w'^). 
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Figure 1. Argument a (solid line) and phase p (dashed line) of 
EJ2Eq depending on (a) distance p and (b) height h. 



From figure 2 can be seen the way in which the 
higher order terms of eq (4) increase and the higher 
order terms of eq (8) decrease with increasing 
distance. The minima of argument in figure la, 
beginning with the chief minimum at about 500 km, 
result from the interference between in — 1)*^ and n^^ 
hop waves and from the interference between the 
first and second order mode, respectively, as can be 
seen from figure 2. This is in accordance with an 
idea from H. Poeverlein [1959]. 

5. Discussion 

The strong equivalence between ray- and mode- 
theory as shown above permits the calculation of 
field-strength values over a wide range of frequency, 
height and distance using eq (4) and eq (8) alter- 
nativeh^ In the case of medium distances it is also 
possible to replace the mode picture by the ray 
picture wlien taking into account higher order 
corrections of the earth ^s curvature and an inhomo- 
geneous and anisotropic ionosphere. In such cases 
mode theory becomes very difficult, and even the 
calculation of the first order mode is extremely 
laborious [Wait and Spies, I960]. In ray theory it 
is sufficient to replace the Sommerfeld ground wave 
by the ground wave above a curved earth according 
to Van der Pol and Bremmer [1949], to exchange the 
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Figure 2. Sum of (a) ray-vectors and (h) mode-vectors depend- 
ing on distance p at constant height h = 70 km. 
The ratio of the scales of figure 2a and figure 2b is 1 : 2. (For details see text.) 



reflection coefficient and to correct for the earth^s 
curvature, for the diffraction [Wait and Murphy, 
1957]. Furthermore, the factor l+/?e(^J must be 
replaced by a more complicated function which is 
valid when 6 is near 90°. 

Examples of calculations covering a wide range of 
frequencies involving the use of both eqs (4) and (8) 
are shown in figure 3 and figure 4. Figure 3 gives the 
argument of Ej2Eo depending on frequency for the 
constant heights A = 70 km (solid line) and A =60 km 
(dashed line). These curves are equivalent to the 
complex transfer ratio of a white noise dipole. 
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Figure 3. Complex transfer ratio of the vertical components of 
the electric field of a white noise dipole for a distance of 1,000 
km. 
Virtual height of the ionosphere A =70 km (solid line) and ^=60 km (dashed 

line). 
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Figure 4. Mean frequency spectrum in the proxiniity of a light- 
ning discharge multiplied by the complex transfer ratio of the 
vertical electric field of a white noise dipole for different 
distances. 

(solid line: h = 70 km; dashed line: /i=60 km) 



These values, iniiltiplied by the mean frequency 
spectrum of a lightning discharge in its proximity 
according to Taylor and Jean [1959], give figure 4a. 
They should be equal to the measured spectrum of a 
mean discharge at 1,000 km distance. Figure 4b 
and figure 4c give similar calculations for distances 
of 2,000 and 3,000 km, respectively. The calculated 
values, especially figures 4a and 4c, show some 
similar features to the mean frequency spectra of 
lightning discharges measured by Obayashi et al. 
[1959], particularly the double maximum between 
10 and 20 kc/s and the relative maximum at about 
40 kc/s. Very similar curves have been presented 
by Wait [1957b] wlio employed the mode sum 
exclusively. 

During a solar (la re the reflection level descends. 
This is represented by a 10 km lower level of h = 60 
km in figure 4 (dashed line). Measurements of 
field strength during such an event generally show 
a decrease below about 15 kc/s and an increase above 
15 kc/s [Obayashi et al., 1959]. In figure 4 tliis is 
only true up to about 30 kc/s. The disagreement 
between calculated and measured valines above 30 
kc/s may be caused by neglecting second orch^r t(M-ms 
of the earth's curvature, which according to Wait and 
wSpies [1900] become more important with increasing 
frequency. 
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